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ABSTRACT 

We present results of deep integral field spectroscopy observations using high resolution opti- 
cal (4150-7200 A) VIMOS VLT spectra, of NGC 4696, the dominant galaxy in the Centaurus 
cluster (Abell 3526). After the Virgo cluster, this is the second nearest (z=0.0104) example 
of a cool core cluster. NGC 4696 is surrounded by a vast, luminous Ha emission line neb- 
ula (Lho = 2.2 X 10^" ergs~^ ). We explore the origin and excitation of the emission-line 
filaments and find their origin consistent with being drawn out, under rising radio bubbles, 
into the intracluster medium as in other similar systems. Contrary to previous observations 
we do not observe evidence for shock excitation of the outer filaments. Our optical spectra are 
consistent with the recent particle heating excitation mechanism of Ferland et al. 

Key words: galaxies: clusters: individual: Centaurus - galaxies: cooling flows - galaxies: 
individual: NGC 4696 



1 INTRODUCTION 

In this paper we report on deep (^-^12 hours) integral field spec- 
troscopy observations of the extended optical emission line system 
surrounding NGC 4696, the brightest cluster galaxy (BCG) in the 
, Centaurus cluster (Abell 3526), which, at a redshift of 2 = 0.0104 
■ is the second nearest example of a 'cool core' cluster. 

Nearly half of all galaxy clusters have dense cores of cool X- 
, ray emitting gas observed as a sharp peak in their X-ray surface 
' brightness profiles. These 'cool core' clusters have short cooling 
times and central temperatures which drop to about a third of the 
cluste r ambient temperature (for a review see IPeterson & FabianI 
I2OO6I and McNamara & Nulsen 200"^. In addition to the sharply 
peaked X-ray surface brightness profiles, many cool core clusters 
have extended filamentary optical emission line systems around 
their central, BCGs (see for examplelHu et al.l 1985l:lHeckman et al.l 
ll989l : ICrawford & FabianI 1992l : lDonahue et al.'l992') and also vast 



only ~ 10 per cent of the predicted star formation is observed. The 
current paradigm is that some form of heating, most likely domi- 
nated by feedback from the central AGN is responsible for regulat- 
ing the heating and cooling in these objects. 

The optical emission line nebulae surrounding the central 
galaxies in these systems, have been found to extend to distances 
of ~I00 kpc in some objects and contain ~10*-10^ M© of 10* K 
gas. The origin and excitation mechanisms of the extended optical 
emission line systems has been much studied but as of yet no one 
mechanism has been found that can explain in detail all the features 
of the emission. 

Mec hanisms proposed incl ude photoionisation by the central 
AGN (e.g. lHeckman et al.ll989r). photoionisation by massive or ex- 
tremely hot stars (e.g. iTerlevich & Melnick 1985; Johnstone et al 
J^8^, heating through shocks (e.g. lCowie et alJll980l ; ISabra et al 



20001) . co nduction or mixing layers in the hot intraclust er medium 



20od : lMcNamara et al.ll9'9^ 



molecular gas and dus t reservoirs (e.g. Edge 2001: Edge et al, 
2OO2I : llaffe et al.l I2OOII : ISalome & Combes! l2003l: iDonahue et"ai 



Without some form of heating, the hot gas in these short cen- 
tral cooling time clusters will have had time to cool and condense 
producing large quantities (10^^-10^'^ Mq) of cold gas and lOOO's 
of solar masses per year of star formation. Although signatures of 
cool and cold gas are observed in the central galaxies, in general 
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(ICM, e.g. lSparks et al.lll989l ; l&awford & Fabi"ai]|l992h. and heat- 
ing of the filaments through suprathermal particles l IFerland et al.l 
[2OO8, 2009). 

The heating and cooling in the Centaurus cluster is very well 
balanced despite the short central cooling time of only 200 Myr. 
NGC 4696 houses a radio source, and multiple bubbles of rela- 
tivistic gas. These are accompanied by soft X-ray filaments and a 
sharp rise in the metal abundance in the centr al 30 kpc, among 
the h i ghest seen in any cluster (~ twice solar, Sanders & FabianI 
l2002ljFabian et al.ll2005l : ISanders & Fabiarj|2006l) . Centaurus also 
has one of the broadest range of measured X-ray temperatures, con- 
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Table 1. Observation log.^ Observations centered on NGC 4696.'^ Scan through galaxy. 
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Figure 1. Our VIMOS [N II]A6583 emission image of NGC 4696 
overlaid with our VIMOS high resolution pointings. The central 
RA and Dec. of box 1 is 12''48'"51^4, -41°18'42.8", box 2 
12''48'"49''.8, -41°18'54.3", box 3 12'»48'"49''.5, -41°18'39.8" and 
box 4 12'' 48"" 48". 0, -41°18'37.3" (all coordinates are in equinox 
J2000). The x and y axes are in pixels. 



taining gas from 0.3 5 to 3.7 keV, over a factor of 10 in temperature 
JSanders et al.'2008^ 

Crawford et al^ ( 120050 presented narrow band Hq-i-[N II] im- 
ages showing the extensive, filamentary Hq nebulosity surrounding 
NGC 4696. This shares the morphology of the soft X-ray filaments 
and of a prominent dust lane. The origin of the complex and ex- 
tensive filamentary system surrounding NGC 4696 has long been 
discussed with suggestions that t he filaments have co oled from the 
ICM, are the result of a merger JSparks et al.ll 19891) or have been 
draw n out into the ICM from the central galax y by rising rad i o bub- 
bles JCrawford et"ai]|2005l) . Recent work bv iFarage et al] J201(]|) 
found in support of the merger origin and presented evidence for 
shock excitation in the filaments. Our data suggest that the fila- 
ments are similar to those seen around many other BCGs many 
of which have not undergone recent merger activity. We uncover 
a second velocity component in the central region which may be 
a filament extending behind the galaxy and discuss the excitation 
mechanisms for the system. 

The observations and data reduction are briefly described in 
§2, analysis, results and a discussion of their implications in §3 
and in §4 we summarise our results and main conclusions. At the 
redshift of the Centaurus cluster {z = 0.0104, 44.3 Mpc) one 
arcsec corresponds to 0.210 kpc (throughout this paper we adopt 
Ho = 71 kms"^ Mpc-^ , ^m = 0.27 and ^a = 0.73). 



2 OBSERVATIONS AND DATA REDUCTION 

Observations were made on 2009 March 27th-30th using the Visi- 
ble MultiC Jbject Spectrograph ( VI MOS) on the VLT in La Paranal, 
Chile (see lLeFevre et alj|2()03l and lZanichelli et alj|2005l for a de- 
scription of the VIMOS IFU and a discussion of data reduction 
techniques). We obtained high resolution orange (HRO), high res- 
olution blue (HRB) and low resolution blue (LRB) data using the 
VIMOS Integral Field Unit (IFU). We used the larger 0.67" fibres 
giving a field of view of 27" x 27" with the HR grism and 54" x 54" 
with the LR grism. Details of the observations are given in Table[T] 
and the HRO pointings are shown in Fig.[T] 

The deep HRO data cover the wavelength of redshifted spec- 
tral lines of [N IJA5199 to [O II] A7330 including five c orona l lines 
the details of which can be found in ICanning et alj ( 120 Uh . The 
HRB data cover the wavelength of redshifted lines of H7 A4341 
to He I A5875. An analysis of the stellar spectra in this object using 
the LRB data will be presented in Canning et al. 201 1 (in prep.). 

The data were reduced by the VIPGu pipeline 
dScodeggio et al.l 120050 . The resulting 3D-datacubes were 
combined with a set of IDL routines (R. Sharp, Private communi- 
cation). 

Instrument flexures in VI MOS are dependent on rotator posi- 
tion and suffer from hysteresis dAmico et alj2008l) . For this reason 
we took calibration frames (three flats and one arc) after each ob- 
servation block. The observation blocks consisted of three science 
exposures. Each were short exposures of only 15 minutes so as to 
minimise extreme flexures. 

The wavelength calibration was done within VIPGI using He, 
Ne and Ar arc frames taken throughout the night. VIPGI performs 
a polynomial fit to the arc lines and we found that the rms devi- 
ations are approximately gaussian with a mean of ~0.05 A and 
dispersion of ~0.02 A. This translates to an error in the velocity of 
~5 kms^^ which is significantly smaller than our FWHM instru- 
mental resolution of 98 km s~^ (a = 42 km s~^ ). 

Flux calibration was done using standard stars observed at the 
beginning and end of every night. Cosmic-ray rejection, final fibre- 
to-fibre transmission corrections, sky subtraction, extinction cor- 
rections and shifting to the object rest frame were performed out- 
side of VIPGI using IDL routines. Sky-subtraction, where done, 
was complicated by the lack of suitable sky fibres and is described 
below. 

Telluric absorption feature corrections for the O2 and H2O 
absorption in the 6000 — 7000 A regime was determined from 4 
observations of standard stars. Standards were observed at the be- 
ginning and the end of the night. Details of the correction are given 
in Appendix lAl 



^ VIPGI- VIMOS Interactive Pipeline Graphical Interface, obtained from 
|http://cosmos. iasf-milano.inaf.it/pandora/ 
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Figure 3. Example spectra with one and two velocity component fits from two pixels south east of the centre of NGC 4696, the pixel coordinates in (x,y) and 
corresponding 'bin' id is shown in the top left hand corner of the plot. Both these spectra are within the spatial region where we have our deepest data. The top 
panel shows the data (black) and model (single gaussian in red, double gaussian in blue), the middle panel shows the separate gaussian components of the two 
velocity component model and the lower panel the fit residuals as data— model; red residuals are data— one velocity component, blue are data— two velocity 
components model. 
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Figure 4. The left hand image shows [N II]A6583 emission. The bright ridge of emission which curls round the south-west of the nucleus is clearly visible. The 
ridge is binned in 2.5 arcsecond squared regions and labelled A to L with A starting near the nucleus. Overlaid are contours from a unshaipmasked broadband 
HST F435w image (shown on the right, both images have the same scale) tracing the dust lanes. The HST image has been smoothed by 20 pixels and 5 pixels 
and the two resulting images subtracted. 



VIPGFs sky subtraction technique is optimised for deep sur- 
vey observations where the field is devoid of extended objects. 
Given this we elect not to use VIPGI's sky subtraction technique 
and instead use specific fibres shown to be lacking in emission 
lines associated with NGC 4696. We median combine the spectra 
to find an average sky spectrum, remove the continuum by means 
of a spline fit and scale and subtract the sky spectra from our object 
fibres. The line profiles are slightly different in each quadrant so 
this process is performed on a quadrant by quadrant basis. 

After basic reduction, transmission correction and corrections 
for sky absorption and emission we bin the spectra in order to reach 
a chosen signal to noise. We try two binning methods: the first bins 
on the basis of the surface brightness of emission in a 2D image by 



following contours of a smoothed image l lSandersll2006f) while the 
second uses Voronoi tesse llations to provide compact spatial bins 
JCappellari & CoDinl2003r) . The results of the two binning methods 
applied to a 2D image of [N II]A6583 emission in NGC 4696 is 
shown in Fig. lA2l 



When binning spectroscopic data one has to be careful not to 
bin over regions which are too extended. The spectral properties are 
unlikely to be similar in two distinct regions of a galaxy and thus 
binning these regions togethe r would lead to spurious re sults. The 
Voronoi tessellation method of JCappellari & CoDinl l l2003h provides 
more compact bins, however the extended optical filament system 
of BCGs adds a further complication. These structures are extended 
in only one direction and as such the requirement for 'round' bins 
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Figure 2. An example BC03 SSP model fit (red) to the observed spectrum 
(black) from bin 5 (60,46), between the wavelength of 4000 - 5200 A. The 
fit residuals are shown in the bottom panel. Any regions where emission 
lines are expected are masked out in the fit. In the above plot the masked 
regions correspond to the wavelengths of, from left to right, H7 A4341, [O 
III]A4363, He II A4686, H/3 A4861, [O III]A4958, [O IIIJA5007 and [N 
I]A5199. 



does not trace their morphology particularly well. Due to the fila- 
mentary structure of the optical line nebulosity surrounding NGC 
4696 w e choose t o bin our spectra using the contour binning tech- 
nique of lSanderj J2006l) . 

Our HRO and HRB spectra overlap between 5000 - 6000 A. 
In order to incorporate all our data we split the cubes into three 
wavelength regions. The regions chosen are 4000 — 5100 A, 
5100 - 5800 A and 5800 - 7500 A. The two HR grisms used 
have a slightly different spectral resolution (see Table [TJ with the 
HRO being marginally larger than that of the HRB grism. In order 
to combine both data sets it is necessary to re-bin the HRO spectra, 
details of the rebinning technique can be found in Appendix IaI 

The strong emission lines of Ha and the [N II] and [S 
II] doublets are fit directly, on a pixel-to-pixel basis across 
the field of vie\y We fit these in IDL using MPFIT ( iMord 
119781 : iMarkwardtl 120091) . The weaker lines are first continuum 
subtracte d using template simple s tellar population (SSP) mod- 
els from iBruzual & Chariot! l l2003h ( hereafter BC03), fit using 
the S TARLIGHT software package dCid Femandes et al.l uOOSl . 
l2009r) . We use basefiles made up of 240 BC03 models spanning 
6 metallicities, covering the range 0.005 — 2.5 Zq, and 40 ages, 
covering the range — 20 Gyr. An example of a fit to the stellar 
continuum between 4000 A and 5200 A is shown in Fig. [2] Af- 
ter continuum subtraction the emission lines are fit using the same 
technique as the stronger lines described below. Where appropriate 
we impose the kinematics determined from the strongest emission 
features. 

The five strongest emission features ([N II]A6548, HqA6563, 
[N II]A6583, [S II]A6717 and [S II]A6730) are fit simultaneously. 
The redshift and velocity dispersion are constrained to be the same 
for one velocity component and the integrated flux of the [N II] 
doublet is tied, the scaling being dictated by the atomic parameters 
dOsterbrock & Ferland 20061). The spectrum is fit between 6590 A 
and 6820 A. A continuum estimate is taken as the region between 
6700—6750 A and subtracted from the spectrum. The continuum 
is fit locally for the Ha+[N II] emission and the [S II] doublet. 




Narrow component 
Broad component - 



, — • — Blueshifted component 
• — Redshifted component - 



Figure 6. The intrinsic Ay extinction in the two velocity components de- 
rived from their Ho/H/3 ratio. The top panel shows the extinction were the 
components to be distinguished as a 'broader' and 'narrower' component, 
the bottom panel as a 'bluer' and 'redder' component. This figure shows 
splitting the two components based on the width of the line as opposed to 
the line centroid is equivalent to splitting the components based on the ex- 
tinction of the emission. 



An example of the fitting technique is shown in Fig. IA5I We fit 
gaussians to the spectra with 1, 2, and 3 velocity components both 
with and without a broad Ha component. We perform an F-test to 
determine which fit to use, this is discussed further in section[3] 



3 RESULTS AND DISCUSSION 

The line emission in the central regions of NGC 4696 is complex 
with many spectra requiring two velocity components to adequately 
fit the observed spectrum. We use the F-test statistic to determine 
which fit to use in the spectral modelling of a particular fibre. The 
F-test allows us to investigate the probability that the data follows 
the simpler of two proposed models nested within each other. We 
use MPFTEST in IDL to calculate the significance of the addition 
of extra parameters in the fit. We choose to accept the additional 
parameters when the probability falls below 0.01. The fits were ex- 
amined by eye to ensure they were reasonable. Fig.[3]shows exam- 
ple fits to two bins; in the first the F-test indicates a single gaussian 
fit is sufficient to explain the data while in the second an additional 
velocity component is required. 

It was not found necessary to fit a broad Ha component. There 
is good agreement between the components found necessary using 
different binning techniques in all areas where the data is of a suffi- 
cient signal to noise to attain a good fit. We adopt, for the remainder 
of the paper, the contour binning method of Sanders (2006) and the 
appropriate F-test results (see Fig. lA4l l. 

The central region of NGC 4696 is exceedingly complex with 
many individual filaments, each of which are likely formed by nar- 
rower interwoven threads as in the case of the extended einissio n 
line system of NGC 1275 in the Perseus cluster ( JFabian et al.|2008h . 
The two velocity components detected in this central region are 
probably due to distinct filaments seen in projection and as such it 
is not obvious that splitting the emission into simple 'blueshifted' 
and 'redshifted' populations will lend any physical insight into this 
system. 

In order to provide a more physical basis for the division of 
the two populations we investigate the emission properties along 
the brightest inner filament. Fig. |4] shows on the same scale the [N 
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Figure 5. To establish a physical basis for distinguishing the two velocity components detected we investigate the kinematic properties were the two compo- 
nents of the gas to be separated by redshift, velocity dispersion or intensity of emission along the bright inner filament. The figures show the most blueshifted 
(top left panel), most redshifted (top right), broadest (middle left), narrowest (middle right), strongest (bottom left) and weakest (bottom right) components. 
Left: The line of sight velocity in kms~^ of each component with respect to the velocity of NGC 4696 (3045 kms~^ ). Right: The FWHM velocity width 
in kms~^ of each velocity component. The (0,0) position corresponds to the nucleus at aRA and Dec. 12''48™49^.28, — 41°18'39.4". 



II]A6583 emission and dust lanes. There is a bright ridge of emis- 
sion curling round the south-west of the nucleus which, in projec- 
tion, coincides exactly with a large dust lane. We bin this region in 
2.5 X 2.5 arcsecond bins (labelled from A to L in Fig.|4]l and investi- 
gate the kinematic properties of the two velocity components were 
they to be seperated based on redshift, velocity width and intensity; 
the results are presented in Fig. [5] 

We would expect, if these two components followed emission 
from two distinct filaments, that the line of sight velocities and the 
FWHM velocity widths would vary smoothly with distance along 
the filament. By which we mean that the velocities are unlikely 
to jump dramatically from bin to bin along the filament length. 
From this argument and using the kinematical data we rule out a 
division based on intensity (the bottom panel of Fig.[5](- To investi- 
gate further separating the components based on redshift and veloc- 
ity width we derive the intrinsic extinction from the Ha/H/3 ratio 
of each component. The results can be seen in Fig. [6] and imply 
that splitting the two velocity components based on the breadth of 
the emission is equivalent to distinguishing them based on intrinsic 
extinction, with the narrower component being more extinguished 
than the broader component in all regions where the Ha and H/3 
lines could be detected at greater than 4 sigma. Splitting the com- 
ponents based on redshift provides no such distinction. From these 
results we suggest that distinguishing the two components by ve- 
locity width as opposed to redshift is the more physical approach. 
The broad component lies closer to the nucleus, coincident with the 
dust lane, whereas the narrow component lies behind, and is there- 
fore more reddened. For the remainder of the paper, unless other- 



wise stated, the two velocity components will be split into 'broader' 
and 'narrower' width components. 



3.1 Reddening 

We dered den the spectra for Galactic extinction using the redden- 
ing law of lCalzetti et alj |2000|) and an E(B-V) value for the stellar 
continuum of 0.11 {Ay ~ 0.34). An extinction to reddening ratio, 
R — Av/E{B — V), of 3.1 is assumed for our Galaxy. This is 
done before shifting the spectra to their rest wavelength. A value 
for the intrinsic reddening is also estimated, where possible, from 
the Ha and H/3 line ratios (Fig. lBlb . 

The logarithmic H/3 extinction coefficient, c, is calculated 
from the observed Ix/Ih^, and intrinsic line ratios /ao//h/3o, 



Ih0 



>-0 ing[/(A)-/(g/?)] 



'H/30 



-10' 



(1) 



assuming a reddening law /(A). The extinction coefficient and red- 
dening can therefore be calculated through, 



A. 



-2.5c/(A). 



(2) 



For the weaker Balmer line we impose the kinematics of the strong 
emission lines and subtract the continuum using template SSP mod- 
els as described in section |2] The intrinsic extinction is calculated 
using the rest wavelengths of the lines and we assume a case B 
intrinsic recombination ratio. This method makes the assumptions 
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Figure 7. [N II]A6583 emission (grey scale) overlaid w itli contours of Ha 
from tlie nan'ow band imaging of iCrawford et alJ bOQSn. The moipliologies 
of other emission lines follow generally that of the [N II]A6583 emission. 



that the dust is an obscuring screen, is homogeneous and that all the 
Ha and H/? emission is from hydrogen recombination, none from 
coUisional excitation. 

The reddening map for a one component velocity fit, on a per 
pixel basis, in the central regions of NGC 4696 is shown in Fig. lBll 
A single component fit is presented as H/3 is weak and the second 
velocity component does not provide a good fit unless the fibres are 
rather coarsely binned, as in Fig. |6] Contours of E(B-V) derived 
from H ST B and I band ima ges and Chandra X-ray Nh column 
density JCrawford et al.ll2005n are overlaid. 

The regions where the value of intrinsic Av is highest corre- 
spond well to the regions where E(B-V) and Nh are greatest though 
the exact position of the peaks are so metimes offset from tho se in 
the X-ray and broad band optical map. ICrawford et alj ( 120051) find 
the regions of greatest extinction have a val ue of Ay ~0.4 which 
is in good agreem ent with previous studies Uorgensen et al.ll 19831 ; 
ISparks et alj 19891) . Our results, f or the same knots, assuming a case 
B intrinsic Ha to H/? line ratio JOsterbrock & Ferlandll2006l) are 
larger than this with A^/'s in the range 1.5-2 for most of the knots 
along the dust lane. In order to make our results consistent with 
these previous imaging studies we require an intrinsic Ha to H/3 
ratio of ~4 (nearly 1 .5 x the case B ratio) requiring a large compo- 
nent of collisional excitation to the hydrogen lines. 

For Seyfert galaxies, radio galaxies and LINERs the in- 
trinsic Ha/H/3 ratio is usually measured to be larger than case 
B, ~3.1, due to the greater importance of collisional excitation 
in regions where there is a harder photoionising spectrum (e.g. 
IVeilleux & Osterbrocklll987D . High ratios have also been seen in 
the ca se of the extended emi ssion line regions surrounding some 
BCGs. lKent & Sargenllf 19791) found a ratio of 4.77 in the filaments 
of NGC 1275, the brig htest cluster galaxy in the Perseus cluster. 
Gemini spectroscopy bv lHatch et al.l ( 120061) confir ms this high ratio 
findin g values of Ha/H/3 of 3-7 in the filaments. IVoit & Donahue! 
( 119971) also find significant reddening in the filaments of Abell 2597 
at 1.5 times the case B ratio. In both NGC 1275 and Abell 2597 the 
authors find no evidence for a deviation from the Galactic extinc- 
tion laws. This implies the properties of the dust in the filaments 
are s imilar to that of our own Galaxy. 

ISparks et al.l jl989r) found, from broad band V and R images, 
that the dust extinction with wavelength behaves 'normally' in 
NGC 4696, that is there is no evidence for a departure from the 
Galactic extinction laws. This makes formation of the filaments 




Figure 8. Soft X-ray emission (0.5-1.5 keV) of the centre of NGC 4696. 
Contours of [N n]A6583 emission are overlaid in black and 5 GHz radio 
contours in white. 



from the hot gas unlikely, and the authors favour a merger origin 
for the dust lane and filaments. More recently it has been suggested 
that the dust lane could be the result of previous dusty star forma- 
tion in the central galaxy which has been drawn out with the cool 
gas into filaments by radio bubbles rising in the ICM. 

To make a co mparison with the extinction corrections of 
iFarage et al.l ( 120101) we examine the extinction in a 3 arcsec- 
ond diameter region approximately 3" north-west of the nucleus 
(taken to be the centre of radio emission and of the X-ray point 
source 12''48'"49\28, -41° 18'39.4"), where we find the peak in 
our Ha emission map. The coordinates of our chosen region are 
(12''48"'49M, -41°18'37.8"). Fitting t he Ha and H , 3 emi ssion 
with a single velocity component, as in iFarage et alj ( l2010n . we 
find the intrinsic Av in this region is 0.33. Using the same values 
of the intrinsic emitted flux ratio (Ha/H/3 = 3.1) as these authors 
gives a total Ay of 0.67. Using the case B ra t io the total Av is 0.91. 
Our value differs from that of IFarage et alj j20ld) by 0.07; this is 
well within the error as a difference in the position of our region of 
1.5" translates to a factor of almost 2 in the Ha/H/3 ratio. 

As the reddening correction varies quite significantly across 
the field of view, due to the presence of the large dust lane, we elect 
to use the local value of the Ha/H/3 ratio where possible to correct 
our fluxes. Where this is not possible we use the median value of 
Av =0.9. 



3.2 Filament morphology 

The detected line emission in our data follo ws closely the emission 
detected in the narrow band Ha imaging of lCrawford et alj ( 120051) 
(Fig. |7j. A central filament is seen with two velocity components 
extending 12 arcseconds to the south-east of the nucleus. A further 
complex system of filaments extending over 25 arcseconds to the 
north-west, south and east are also detected. Only the north-east 
appears devoid of filaments. All emission lines, with the exception 
of [O III], where the lines were strong enough not to impose the 
kinematics of the Ha emission were found to have the same mor- 
phology as the Ha line suggesting the emission occurs from the 
same gas. There also appears to be dust lanes present in all fila- 
ments except the weak southern-most filament (Fig.|4ll. 
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Figure 9. The morphologies of the broader (left) and narrower (right) velocity components. The peak in the emission between the two components is offset 
by ~2". The broader component has a much more circular profile while the narrow component has a steeper central peak and exhibits a less smooth structure. 
The two images have been smoothed by a gaussian with FWHM 3.1" (cr =1.3") and are plotted on the same spatial scale; the cyan 'X' marks the nucleus. The 
intensity scale has been adjusted individually to achieve the best contrast. 




Figure 10. The [O III]A5007 (left) and [N II]A6583 (right) emission in the central region of the galaxy. The peak in emission is offset by ~4" between the 
two images. The two images have been smoothed by a gaussian with FWHM 3.1" {a =1.3") and are plotted on the same spatial scale. The intensity scale has 
been adjusted individually to achieve the best contrast. 



The emission line gas traces well the inner regions of the soft 
X-ray filaments (Fig.[8](. However, the soft X-ray filaments extend 
much further up to the north-east of the galaxy leading towards 
a depression in the X-ray gas thought to be a ghost bubble, of 

S ilasma from the radi o source, risen bouyantly from the nucleus 
Crawford et al.ll2005r) . The area just north-east of the nucleus, de- 
void of both optical emission line gas and soft X-ray filaments co- 
incides with a prominant radio bubble. 

The two velocity components have a different morphological 
structure. Fig.|9]shows the [N IIJA6583 line emission in this central 
region, smoothed by a gaussian with FWHM 3.1" (a =1.3", 2 pix- 
els). The left hand image is the emission from the broader compo- 
nent and shows two clear peaks in emission. The weaker peak to the 
east coincides with the centre of radio emission and the weak X-ra y 
point source (12'M8'"4r.28, -41° 18'39.4". iTavlor et ai]|2002h : 
throughout this paper we define the nucleus as this point. 5 arc- 
seconds to the west is the point of largest flux for the broad emis- 
sion component. The surface brightness contours show the emis- 
sion about this point is smooth and circular in the central 5 arcsec- 
onds before stretching towards the south-east of the nucleus. The 
right hand image shows the narrow emission component which is 
far less smooth in morphology. The emission peaks in a bar stretch- 
ing north-west from the nucleus then curls clockwise round the nu- 
cleus. We find the contours of surface brightness of both of these 
components deviates from that of the dust near to the nucleus as 
found bylfarage et al. (20101). 

Fig.[TO]shows a single component fit to the [O IIIJA5007 emis- 



sion (left), beside an [N IIJA6583 single component fit emission 
map (right). The [O III] emission is weakly detected in the inner 
filament (we do not have a HRB pointing in the outer regions) but 
very strong near to the centre of radio emission. The second peak 
of [O III] emission to the west which coicides with the peak in the 
other emission lines is much weaker, implying much of the [O III] 
emi ssion in the central re gions is associated with th e AGN. 

ISparks et al.l l ll989l) and iFarage et al.l 1 I2OIOI) suggest NGC 
4696 has an intrinsically one sided spiral structure due to the ap- 
parent close association of the main dust lane and line emission. 
While we also observe a very striking close conespondence be- 
tween the optical filaments and the many dust lanes of NGC 4696 
(see Fig.|4ll, our detection of at least two velocity components to the 
optical emission line gas, exhibiting different extinction properties, 
provides evidence contrary to this suggestion. NGC 4696 is clearly 
a very complex galaxy surrounded by an intricate filamentary sys- 
tem. The high extinction of the narrow velocity component and its 
less smooth morphology is interpreted as a filament located behind 
the main, more centrally situated broader emission line component. 



3.3 Emission line fluxes 

The spectrum of a region with a radius of 10 arcsec centred on the 
radio core, with and without continuum subtraction, is shown in 
Fig-El with the strongest lines marked. The spectrum is similar to 
other filament emission in cool core BCGs, showing strong low ion- 
isation emission lines. NGC 4696 has stronger [O III]A5007 emis- 
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Figure 11. Top: The spectnjm of NGC 4696 within a radius of 10 arcseconds, sky lines that are confused with object emission are not subtracted. Bottom: The 
same spectrum after continuum subtraction using BC03 models. Object emission lines and confused sky lines are labelled. 



Line 






Total flux 
ergcm"^ s^^ 






Total luminosity 
ergs-i 




Flux / Fhc 


H7 A4341 






< 2.2 X 10-15 






< 4.9 X 10^8 




0.02 


[OIII] A4363 






< 1.5 X 10-15 






< 3.3 X lO^s 




0.01 


He II A4686 






< 2.4 X 10-15 






< 5.5 X 10^8 




0.02 


H^ A4861 


2.0 


X 


10-1'' ±9.9 X 10" 


-16 


4.5 


X 10^3 ± 2.3 X 


1038 


0.20 


[OIII] A4958 


9.2 


X 


10-15 ±1.2 X 10- 


-15 


2.1 


X 1039 ± 2.8 X 


1038 


0.09 


[OIII] A5007 


2.8 


X 


10-" ±5.5 X 10- 


-IB 


6.3 


X 1039 ± 1.3 X 


1038 


0.28 


[NI] A5199 


5.4 


X 


10-15 ±1.3 X 10" 


-15 


1.2 


X 10^9 ± 2.9 X 


1038 


0.06 


[CV] A5309 






<6.0x 10-16 






< 1.4 X 10^8 




0.01 


[OI] A5577 






< 2.9 X 10-15 






< 6.6 X 1038 




0.03 


[Nil] A5755 






< 4.2 X 10-15 






< 9.6 X 1038 




0.04 


He I A5875** 


1.6 


X 


10-" ±1.4 X 10- 


-15 


3.8 


X 10^9 ± 3.3 X 


1038 


0.17 


[OI] A6300** 


1.7 


X 


lO-" ±4.0 X 10- 


-15 


4.0 


X 10^9 ± 9.0 X 


1038 


0.18 


[OI] A6363** 


5.7 


X 


10-15 ±3.0 X 10- 


-15 


1.3 


X 1039 ± 7.0 X 


1038 


0.06 


[Nil] A6548 


5.4 


X 


10-" ±1.3 X 10- 


-14 


1.2 


X 10*" ± 3.0 X 


1039 


0.56 


Ha A6563** 


9.8 


X 


10-" ±5.9 X 10- 


-15 


2.2 


X 10*0 ± 1.4 X 


1039 


1.00 


[Nil] A6583 


1.8 


X 


10-1-' ±9.5 X 10- 


-15 


4.2 


X 10"" ± 2.2 X 


1039 


1.87 


[SII] A6717 


5.5 


X 


10-"±5.7x 10- 


-15 


1.3 


X 10« ± 1.3 X 


1039 


0.56 


[SII] A6731 


4.5 


X 


10-" ±5.2 X 10" 


-15 


1.0 


X 10*'' ± 1.2 X 


1039 


0.46 


[Call] A7291** 






< 1.7 X 10-" 






< 3.9 X 1039 




0.18 


[Oil] A7320** 






< 3.7 X 10-" 






< 8.4 X 1039 




0.38 


[Call] A7325** 






<3.7x 10-" 






< 8.5 X 1039 




0.38 


[Oil] A7330** 






< 2.1 X 10-" 






< 4.8 X 1039 




0.21 



Table 2. Total fluxes and luminosities of emission lines in NGC 4696 after correction for extinction using the median Ay = 0.9. Emission lines blueward of 
5100A are only measured in the inner 27" as we do not have HRB data beyond this. Where the lines are marked by < the value indicates a 90 per cent upper 
limit (Ax^ = 2.7). Where the Unes are marked with ** the object emission was confused with other emission or sky Unes and it was necessary to fit these 
simultaneously. All values are determined after subtraction of the continuum. 



sion than found in NGC 1275 or Abell 2597 though the majority 
of this emission coincides with the nucleus and is likely associated 
with the AGN. 

We compare the dereddened flux in Hq emission to that of 
ICrawford et al.l ( 120051) . Using narrow band Ha imaging the authors 
find a surface brightness in the inner (^ 8" from the radio core) 
filaments of 3.5 x 10-^^^ ergcm-^ s-iarcsec-^ and in the fainter 
outer filaments (~ 27") of 1.4 x 10^^® ergcm-^ s-iarcsec-^. 
They obtain a lower limit on the total Hq luminosity beyond a ra- 
dius of 3.5" of 1.5 X 10*° ergs-i . Our data are consistent with 
this with the average surface brightness of the filaments 8" and 27" 



from the radio core being 2.2 x IQ-i® erg cm-'^ s~'^arcsec-^ and 
1.3 x 10^^'' ergcm-^ s-iarcsec-'^ respectively. We obtain a to- 
tal extinction corrected Hq luminosity beyond a radius of 3.5" of 
1.7 x 10*" erg s~^ and a total extinction corrected Hq luminosity 
of 2.2 x 10*° erg s~^ . Our derived Hq and [N II] fluxes and lumi- 
no sities are slightly la rger than, though still consistent with, those 
of iFarage et al.l ( 120101) . This is expected as our HRO data probes 
more of the filament system. 

The total luminosities of all detected lines and the 90 per cent 
upper limits of emission lines not found in our spectra are given in 
Table |2] Emission lines are fit as described in Section |2] although 
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Figure 12. Regions A and B, both 3 pixels in radius. Region A is coincident 
with the nucleu s defined as the cen tre of radio emission (12'' 48™ 49". 28, 
-41°18'39.4" )Tavlor et al] |20QA. while region B is --12 arcseconds to 
the south-east of the nucleus, on the outer edge of the inner filament. 
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Figure 13. The non-continuum subtracted HRB spectra in regions A (top 
panel) and B (bottom panel). The [O in]A5007/H/3 A4861 ratio varies dra- 
matically from ~2 (region A) to '^0.7 (region B, see table[3). 



some lines are confused with other emission lines or sky and there- 
fore warrant special attention. He I, [O I]A6300, Ha, [Ca II] and 
emission lines of [O 11] suffer from confusion with sky emission, 
in each of these cases the sky emission is fit simultaneously with 
the object emission lines. The red wing of the object [O I]A6363 
e mission coincides with the [Fe X]A6374 emission line discussed 
in lCanning et alj ( 1201 ih . The flux of this line is linked by atomic 
physics to the flux of the [O I] A6300 emission and we fit these, the 
sky [O I] lines and the [Fe X] A6374 emission line together. We only 
have HRB data from the central ~27" of NGC 4696, therefore all 
emission lines blueward of 5100A may be lower limits. 



3.3.1 [O HI] emission 

A characteristic of the extended optical emission filaments in many 
BCGs is their low [O III] emission with respect to the hydrogen 
recombination lines (see for example iHatch et al.ll2006l) and sur- 



Line 


F\/FuaA 


Fa/ Fhc b 


H/3 A4861 


0.37 


0.23 


[OIII] A4958 


0.26 


- 


[OIII] A5007 


0.77 


0.15 


[NI] A5199 


0.23 


0.19 


He I A5875** 


0.15 


- 


[OI] A6300** 


0.18 


0.21 


[OI] A6363** 


0.06 


0.07 


[Nil] A6548 


1.04 


0.57 


Ha A6563** 


1.00 


1.00 


[Nil] A6583 


3.10 


1.70 


[SII] A6717 


0.92 


0.39 


[SII] A6731 


0.58 


0.33 



Table 3. Extinction corrected F^/ Fna in regions A and B (see Fig. I12t . 
A coincides with the nucleus, defined as the centre of radio emission while 
B is a point ~12 arcseconds away on the north-eastern edge of the inner 
filament. Where the lines are marked with **, the object emission was con- 
fused with other emission or sky lines and it was necessary to fit these si- 
multaneously. The measured Ha fluxes in the two regions are: Fjja A = 



5.29x10 ^^ erg cm ^s ^andFno 



: 2.99 xlQ- 



erg cm 



prisingly high [Ne III] emissio n lines (see discussion of charge ex- 
change in lFerland et al.ll2009f) . The ionisation potential of Ne^ is 
greater than that of 0+ so if we observe [Ne III] lines we would 
expect to see those of [O III] as well. 

The total [O III]A5007 emission in NGC 4696 is fairly strong 
with an [O III]A5007/H/3 A4861 intensity ratio measured as 1.4 in 
the inner 27x27 arcseconds (the field-of-view of our HRB grism). 
However, as discussed in Section [T2l the surface brightness of [O 
III] emission is steeply peaked and located near to the centre of 
radio emission. It is therefore probable that much of the [O III] 
emission is due to the central AGN. 

We investigate this possibility by choosing two regions; one 
coincident with the radio nucleus (region A) and one offset by 
~12 arcseconds, on the farthest edge of the inner filament from 
the nucleus (region B). Spectra are extracted in both regions from 
a 3 pixel (2") radius circle, the regions are shown in Fig. [12] 
and the spectral regions containing the H/3 line and [O III] emis- 
sion are shown in Fig. [13] The ratio of [O ni]A5007/H/3 A4861 
clearly varies in the two regions with the central region having 
a ratio of ~2 and the filament emission having a much lower [O 
III] A5007/H/3 A4861 ratio of only ~0.7. Ratios of the fluxes of the 
detected lines to Ha for the two regions are shown in Table|3] 

Observation s of t he inner r egions of NGC 1275 
JKent & Sargent 1 19791 : ISabra et al.l I2OO0I) have found [O 
III]A5007/H/3 A4861 ratios of 0.6 in the inner 24 arcsec (8 kpc), 
however the [O III]A5007 emission in the outer filam ents, beyond 
10 kp c, is not detected in the long slit spectra of JHatch et al.l 
I I2OO6I) . The authors are able to put an upper limit on the [O 
III] A5007/Ha A6563 emission in a region 18 kpc from the nucleus 
of <0.03, their limit on [O III]A5007/H/3 A4861 is <0.13. 

Using long sl i t spec tra from the INT and FOS spectra 
jjohnstone & FabiatJ 119881) found that the [O III]A5007/H/? A4861 
ratio decreases away from the core in both NGC 1275 and Abell 
1795. In NGC 1275 their spectra extend to a radius of 8 kpc where 
they find the ratio drops by a third from the central value to about 
0.5. Similarly in Abell 1795 the ratio peaks at about 1 in the nu- 
cleus and decreases to ~0.6 at the outer extent of the slit which 
correspond to a radius of 6 kpc. Low [O III]A5007/H/3 A4861 ratios 
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Figure 14. Top: (Left) The line of sight velocity relative to the velocity of the galaxy, 3045 kms~^ (z=0. 010157), from lPostman & Lauen 1199511 . if a single 
component fit is made over the whole field of view and (right) the line of sight velocity shown only where one velocity component was found to be sufficient. 
The figures shown are fit from the strong [N II]A6583 emission line. Middle: The same as the above for the two component fit. The left hand plot indicates 
the broader component and the right hand plot the narrower component. Bottom: The blueshifted (left) and redshifted (right) velocity components. The colour 
bar is in units of km s~^ and the axis (0,0) position con'esponds to RA and Dec of 12''48'"49'' .28, — 41° 18'39.4", coinciding with the core of the radio 
emission. 



have also been observed in th e vast emiss ion line nebula associated 
with the BCG of Abell 2597 dVoit & Donahue 1997) . 

Shock models tend to predict higher [O III]A5007/H/3 A4861 
intensity ratios and larger fluxes of [O III]A4363 than those seen 
in the o uter extents of the emission line nebulae s urrounding these 
objects JSabra et al.lF2000l ; IVoit & Donahue! 1 19971) . We obtain only 
upper limits on the [O III]A4363 emission across the field-of-view 
with the upper limit on the [O IIIJA4363/H/3 A4861 intensity ratios 
of ~0.08. Unfortunately we only acquired one pointing of HRB 
data in NGC 4696 so are unable to investigate the ratio in the very 
outer extent of the nebulae. 

Our spectra do not extend far enoug h into the blue to ob serve 
the [Ne III]A3869 optical line, however, iFarage et alj ( 1201(31) have 
detected this at the 3(t level (see their Fig. 3). The infra-red line 
of [Ne III]A15.56 ^tm has also been clear ly detected using Spitzer 
IRS SH spectra bv ljohnstone et al.l l l2007h . Photoionisation models 



for the excitation of the emission line nebulae struggle to account 
for the anomalously low [O III] yet high [Ne III] emissio n unless 
coupled with high metallicities. I Sanders & FabianI ( l200q) find the 
oxygen and Ne abundances to peak at about solar towards the centre 
of the Centaurus cluster. However, the measured Ne abundance is 
highly dependent on the spectral model and varies between 0.5-2 
between the SPEX and APEC models. 



^Particle ionisation models put forward bv lFerland et al.l 1 120081 

120090 have so far been successful at reproducing the observed ra- 
tios due to the importance of charge-transfer. Using CLOUDY sim- 
ulations (last described bv lFerland et al.lll998i) the authors find the 
charge-transfer recombination of O"*" and O"*"^ is very fast com- 
pared to He^, Ne^ and Ne^'^ resulting in significant [Ne IT] and 
[Ne III] emission compared with [O III]. 
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Figure 15. Top: (Left) The velocity width at full width half maximum if a single velocity component fit is assumed over the whole field of view and (right) the 
velocity width where only one velocity component was found to be sufficient. The figures shown are fit from the strong [N II]A6583 emission line. Middle: 
The same as the above for the two component fit. The left hand plot indicates the broader component and the right hand plot the narrower component. Bottom: 
The blueshifted (left) and redshifted (right) velocity components. The color bar is in units of km s~^ and has been corrected for the instrumental width. The 
axis (0,0) position coiTesponds to RA and Dec of 12'' 48'" 49^^.28, — 41°18'39.4", coinciding with the core of the radio emission. 



3.4 Velocity structure 



Maps of the kinematics of the optical emission line nebulae are 
shown in Figs. [14] and [15] The line-of-sight velocities are given 
with respect to the redshift of NGC 4696 dete rmined to be 2 = 
0.010157 (3045 kms"^) bv Postman & Lau er (1995). The ve- 
locity dispersions are presented as the FWHM velocities measured 
from the gaussian fits to the strongest emission lines in our spec- 
tra and have been corrected for the measured instrumental width, 
determined from fitting sky lines close to the observed Ha emis- 
sion line (FWHM 98 kms"^ , a = 42kms"^ ). All emission 
lines strong enough t o map were found to have essentially the same 
velocities, as seen in lFarage et al.l ( I2OI0I) . therefore we only show 
mapsof [NIIJA6583. 



3.5 Central velocity structure 

The line-of-sight velocity of the central region of the galaxy 
changes smoothly by ~400 krns~^ across the north to south di- 
rection with the southern region being more redshifted, perhaps in- 
dicating a component of rotation of the nebula about the nucleus. 
This gradient in velocity structure is seen in both the broad and 
narrow velocity components but to different extents (see Fig.|16t. 

The broad component has a wider distribution of line-of- 
sight velocities, with the northern-most and southern-most veloc- 
ities varying by ^^500 kms^^ while in the narrow component 
the velocity changes by ~250 kms^^ from the northern-most to 
southern-most regions. This limited range in velocities is similar 
to those observed in the extended filaments seen to the north-west, 
south and east of the nucleus. 

The broad velocity component has a FWHM of ^^300— 400 
kms^^ while the narrow component has widths in the range 
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Figure 16. The distribution of line-of-sight velocities in the filaments and 
cental regions. 
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Figure 17. The line-of-sight velocities in the outer filaments. The velocity 
zero-point is taken as 3045 kms"'^ . 



50—200 kms^ . The line widths of both components are much 
greater than the thermal widths of the hydrogen gas which is broad- 
ened by ~10 km s"'^ at lO"* K implying either the gas is broadened 
by shocks or by turbulence. 



3.6 Outer filaments 

The distribution of line-of-sight velocities of the outer filaments 
is less than the central region, with a range of ~200 kms~^ in 
most cases. We bin the outer emission in 3 by 3 arcsecond bins 
to examine the velocity structure more closely along the length of 
the filaments, the results are shown in Fig. [Tt] The line-of-sight 
velocity zero-point is taken as 3045 kms^^ . 

The eastern filament is the dominant structure curling round 
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Figure 18. The FWHM velocity widths in the outer filaments after correc- 
tion for the instrumental width. 



the south-east of the galaxy, coinciding closely wit h the outer dust 
lanes a nd soft X-ray filaments. The Ha imaging of ICrawford et al.l 
( 120051) shows this filament is in fact m ade of several thinn er threads 
such as those resolved in NGC 1275 JFabian et al1l2008l) . The fil- 
ament extends over 30" (6 kpc) from the nucleus but exhibits a 
very narrow distribution in line-of-sight velocities peaking at the 
same velocity as the galaxy (see Fig. |16t . The low velocities sug- 
gest we are vewing this filament close to the plane of the sky. Both 
the eastern filament and the north-west filament appear to curve 
round the radio bubble. The north-west filament has a larger varia- 
tion in line-of-sight velocities, though the variation is still smooth, 
with the innermost part more blueshifted and the outermost part 
more redshifted. 

The southern filament exhibits a divide in line-of-si ght veloc- 
ity alo ng the filaments length such as was found by iHatch et al.l 
( 120061) in the northern filament of NGC 1275. The northern most 
part, closest to the nucleus, displays a redshift of ~140 kins~^ 
whereas the region farthest from the nucleus displays a blueshift 
of ~ —50 kms^^ . The regions closest and farthest from the nu- 
cleus along this filament may thus be moving in opposite directions 
with the filament either collapsing in or stretching out, however the 
filament is clearly curved and this scenario is complicated by pro- 
jection effects. 

The velocity widths of the emission lines in the outer filaments 
are fairly low, ranging from 50—250 kms^^ with the more north- 
ern filaments having a slightly larger velocity width than the south- 
em ones. The range of velocity widths measu red in the se filaments 
is similar to those seen in the Perseus cluster ( IHatch et al..,2006,) . 



3.7 Implications for filament origin 

Kinematically the central narrow emission line component bears 
more similarity to the outer filaments than the broad component 
having both a lower velocity dispersion and a narrower distribu- 
tion of line-of-sight velocities. This supports the interpretation of 
the narrow component being a separate filament farther from the 
nucleus, seen in projection with the main central broad component. 
This would suggest the filamentary system is not an intrinsically 
one-sided spiral structure and it is unclear how a merger origin 
could account for the diversity in morphology and kinematics of 
these filaments. 

If the filaments had an inflow origin, as predicted by early 
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Figure 19. lonisadon plots JBaldwin et"ai]ll98lh of the outer and inner regions in tlie central 27" of NGC 4696. Top: [O III]A5007/H/3 A4861 verses 
[N II]A6583/Ha A6563 emission. Middle: [O III]A5007/H/3 A4861 verses [O I]A 6300/Ha A6563 emis sion. Bottom: [O III]A5007/ H/3 A4861 verses [S 
II](A6716 + A6731)/Ho A6563 emission. Solid lines are from lKewlev et alJJ2006l) and dashed lines from lOsterbrock & Ferlandl ilOOfil) . 



cooling flow models JFabian et al.l Il984h . we would expect that 
there would be a negative radial gradient in the widths of the emis- 
sion lines due to bo th greater acceleration towards the centre and 
line-of-sight effects Meckman et alJll989l) . This is not observed in 
any of the outer filaments with the velocity dispersion being largely 
unvarying with the exception of the north-west filament where the 
velocity dispersion appears erratic along the filament. This filament 
exhibits knots in the Ha image, just north-east of the bend in the fil- 
ament and again at the end. Another fainter filament is also present 
in the Ho image, which crosses or joins the north-west filament at 
this bend where the measured velocity width is largest (see Fig.llSt. 
It is therefore likely that the measured velocity width in this region 
is large and erratic due to two or more filaments overlapping, in 
the line-of-sight, with slightly different velocities and resulting in a 
broa d, strong peak o f emiss ion. 

ICrawford et al.l ( 120050 suggest the filamentary system is the 
consequence of buoyantly rising radio bubbles which draw fila- 
ments of cool and cold gas and dust up be neath them as has also 
been seen in other brightest cluster galaxies JBohringer et al. 



Churazov et"aLl 1200 lE fRevnolds et al] l2005riFabian et al.l 
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Hatch et al.ll2006f) . This is supported by a depression in the X-ray 



intensity and thermal pressure maps to the north-east of the nucleus, 
just beyond the optical and soft X-ray filaments (see their Fig. 5). 
A similar depression is seen in the south-west and low frequency 
(330 MHz) radio emission also extends towards both these struc- 
tures. The kinematics of the filaments are consistent with this pic- 
ture with smoothly varying line-of-sight velocities and low velocity 
dispersions which show little variation across the whole length of 
some filaments. 



3.8 Emission line ratios 

The emission line ratios in the optical nebulosity are important 
tracers of the excita tion processes dominating their behaviour. 
iHeckman et al.l ( 119890 suggested there was a bimodal population of 
BCGs based on their [N II]/Ha ratios and Ha luminos i ty. Ho wever, 
using a much larger cluster sample ICrawford et al.l ( 119991) have 
shown that this is in fact a continuous spectrum with all objects 
sharing the charact eristic low-ion i sation emission line spectrum. 
In a study by IWilman et al.l ( l200q) of 4 high Ha luminosity 
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Figure 20. [N II]A6583/Ha A6563 ratio plotted against the FWHM veloc- 
ity with of each pixel and against the Ha flux. Only pixels where both [N 
II]A6583 and Ha A6563 were detected at the 4 a level are plotted. Colour 
codeing same as in Fig. 1191 



BCGs all situated in cool core clusters and with the same LINER 
emission properties, the authors found no evidence for variations in 
line ratios as a function of either position or Ha luminosity across 
the galaxies. They interpret this as evidence for both a single ioni- 
sation source for the line emission across each galaxy and a single 
domin ant excitation mechanism in all of their galaxies. iHatch et al.l 
( 120071) found, in a similar study of optical nebulosities surrounding 
6 BCGs spanning a range of optical. X-ray and radio properties, 
that the ionisation state of the optical nebulosities were not uni- 
form and so concluded the converse; that a single scenario cannot 
be responsib le for the optical line emission in their sample. A re- 
cent study bv lEdwards et al.l ( 120091) of a further 9 clusters, including 
in their sample both cool and non-cool core clusters, found diverse 
emission line ratio s and morphologies. This supports the finding by 
iHatch et al.l 1 120070 that there is no one clear consistent mechanism 
to explain the origin of the line emission or to separate BCGs in 
cool and non-cool core clusters. Despite this there are some clear 
similarities in that a relatively hard ionising source is needed and 
that the systems are always well described by a LINER-like emis- 
sion line spectrum. 

iHatch et alj ( l200q) studied in detail the spectral features in 
NGC 1275, a similar object in the centre of the Perseus cluster, and 
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Figure 21. [N II]A6583/Hq: A6563 emission ratios in the outer north-west 
and southern filaments. The emission is binned in 3"x3" regions across 
these filaments. Only regions where both emission lines were detected at 
greater than 4 sigma are shown. 



found radial variations in the [N n]/HQ emission line ratio along 
the filaments with the inner regions having higher [N Il]/Ha ra- 
tios. This could be due to an extra energy source closer in causing 
more heating per hydr ogen ionisa t ion or due to metallicity varia- 
tions in the filaments. ISarzi et al.l 1 120061) also found variations in 
ratios of [O Il]/H/3 along the length of the filaments in M87, the 
central galaxy in the Virgo cluster. 

Maps of our em ission line ratio s, binned using the contour 
binning algorithm of ISandersI ( llOOq) are shown in Fig. IB2I The 
top panel showing the single gaussian fits, the middle showing the 
broad central velocity component and the bottom panel in each case 
shows the narrow component. Only pixels where the line emis- 
sion is detected at g reater than Aa are shown. A ionisation plot 
JBaldwin et al.lll98lh of the same pixels is shown in Fig.[T9] Black 
points indicate the pixels where one emission line component was 
found sufficient to model the emission lines, red points correspond 
to the broad lines and blue to the narrow. The [O IIl]]/H/3 ratios 
cannot be traced across the whole length of the outer filaments as 
our HRB data only cover the central 27 arcseconds, as such we 
cannot plot the most extended parts of the outer filaments on Fig. 

m 

The scatter in Fig. [19] is large, however NGC 4696 clearly 
exhibits LINER-like emission in both the central and extended 
regions. The broad inner component has in general a higher [N 
II]/Ha ratio than the narrow or outer regions, which are virtually 
indistinguishable from each other. If the broad component is from 
more c entrally located em ission this is consistent with the observa- 
tion of lHatch et all ( |2006|) in NGC 1275 that [N II]/Hq ratios are in 
general higher towards the centre of the galaxy. If the filaments are 
shock excited, we should see a gradient in emission line ratios with 
velocity dispersion of the gas. Fig.[20]shows this ratio for the same 
three regions as in Fig. [19] No correlation is seen in [N Il]/Ha ratio 
with velocity width with most outer filaments (black points) having 
FWHM velocity widths of 100-200 kms^^ but ratios which vary 
from 1.8-3.3. We now examine the variation in emission line ratios 
across the outer filaments in more detail by binning the emission in 
3" by 3" regions across the field-of-view (see Fig.l2U. 
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Figure 22. Variation in emission line ratios as a function of projected dis- 
tance along the north-west and southern filaments. Regions where spectra 
were extracted are shown in Fig. 12 II 
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Figure 23. [N II]A6583/Ha A6563 emission ratio as a function of FWHM 
velocity width in the north-west and southern filament region. 



3.9 Outer filaments 

The north-west and southern filaments show a difference in the [N 
II]A6583/Ha A6563 emission ratios across their projected length 
with the southern filament decreasing and the north-western fila- 
ment increasing in [N II]A6583/Ha ratio (Fig.l21bnd Fig.l22t. The 
innermost regions of both filaments have a ratio ~2, the average 
ratio of the galaxy. This difference in va riation is in c ontrast to that 
seen in the filaments of NGC 1275 by JH^cheTal] ( I2OO61) where 
there is a clear radial trend in [N II]A6583/Ha intensity ratio seen 
in all long slits, with the ratio decreasing with projected distance 
from the nucleus. 
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A higher [N II]A6583/Ha intensity ratio indicates that there 
is more heating occurring in the gas per hydrogen ionisation. The 
difference in emission line ratios in these two filaments could be 
due to the proximity of the filaments to the expanding radio bubble. 
The arm of the north-western filament appears to encircle the radio 
bubble while the southern filament appears spatially less connected 
with the current radio emission. 

The [N II]A6583/Ha intensity ratio is also sensitive to varia- 
tions in the metallicity of the gas . X-ray observations of the Centau- 
rus cluster JSanders et alJuOOSf) have shown there is a high abun- 
dance of nitrogen emission in the central regions. Oxygen is an ef- 
ficient coolent and high oxygen abundances in the gas would act 
to decrease the gas temperature and so reduce the observed [N 
II]A6583/Ha ratio, however the nitrogen/oxygen ratio in the hot 
gas is also high. The bottom panel of Fig.[22]shows the trend in [N 
II]A6583/[0 I]A6300 emission line ratios in these same filaments. 
There is weak evidence for a decreasing trend in [N II]A6583/[0 
I]A6300 ratio in these filaments. 

The variation of [N II]A6583/Ha emission line ratio with 
FWHM velocity width of the gas, in these two outer filaments, is 
shown in Fig. [23] There is no correlation between the measured 
ratios and the width of the lines in either of the southern or north- 
western filaments. The majority of the emission has velocity widths 
less than 250 kms~^ and a emission line ratio of ~1. 8-2.5. The 
two very broad regions in the north-west filament correspond to 
bright knots of emission which could be formed by the intersection 
of two or more individual filaments which cannot be resolved in 
our data. 



4 CONCLUSIONS 

The origin and excitation of the complex extended emission line 
nebula in BCGs and more specifically in NGC 4696 has been a 
matter of some controversy. Fig. |4] shows the close correspondence 
between the dust absorption and the [N II] emission tracing the 
10* K emission-line nebula in NGC 4696. Virtually all extended 
filame nts are spatially coincident with dust emission. | Sparks et al.| 
( 119891) do not find any deviations in the dust extinction behaviour 
with wavelength with respect to dust within our Galaxy, similar 
results have been found in other BCGs. The authors suggest that 
this favours a merger origin for the emission line filaments over the 
formation of the filament system straight from the hot gas. 

There is however, a third possibility; that the gas originated 
due to dusty star formation in AGB atmospheres in the brightest 
cluster galaxy and has been drawn out, rather tha n falling in, by 
buoya ntly rising radio bubbles inflated by the AGN JCrawford et al.l 
120050 . This scenario is supported by Chandra X-ray observations 
that reveal pressure minima towards the north-east at the end of 
the soft X-ray filaments and another minima towards the south- 
west. These minima coincide with the outer edge of the 330 MHz 
radio emission and app ear to be surrou nded by regions of higher 
metallicity X-ray gas ( iFabian et al.ll2005i) . 

Our deep optical IFU data allow us to explore in detail the 
morphology, kinematics and ionisation properties of the optical 
nebulosity surrounding NGC 4696. Our observations show: 

(i) there are at least two velocity components to the emission 
line gas in the inner regions of the galaxy. We interpret these as 
'broader' and 'narrower' velocity width components, 

(ii) these components have a different morphology, with the nar- 
row component having a less smooth morphology and a striking bar 
of bright emission while the broad component exhibits two peaks 



in emission, one peak coincident with the radio core and the other 
4 arcsec to the west (see Fig.[9l(, 

(iii) the peak in [O III] emission is spatially coincident with the 
nucleus and is offset from the brightest peak in [N II] emission (see 
Fig.m, 

(iv) the kinematics and emission line spectra of the extended 
filaments are similar to filaments surrounding other BCGs, 

(v) we do not find any evidence to suggest there is a correlation 
between the velocity width and [N Il]/Ha ratio in the filaments of 
NGC 4696, 

(vi) there is very little [O III] emission in the outer filaments as 
observed in similar systems. 

The central 12 arcseconds radius in NGC 4696 contains at 
least two components to the emission line gas, we interpret these 
as a broad velocity component and a second narrower component. 
The reddening and kinematical properties of this second compo- 
nent lead us to suggest it is a filament with properties much like 
the other extended filaments, that is located behind the main galaxy 
and thus more obscured by the intervening dust lanes. The mor- 
phology of the central region was previously thought to be a simple 
one-sid ed spiral structure perhap s caused by the tail of an infalling 
galaxy JSparks et a l. 1989: Farage et al.l2010r) however the filament 
system is clearly complex and it is not immediately obvious how a 
merger could account for the diversity in morphology of the fila- 
ments. 

We find smoothly varying line-of-sight velocities in the outer 
filaments and low velocity widths (FWHM 50-250 kms^^ ), with 
little variation along the length of many filaments, similar to those 
found i n NGC 1275, the brightest cluster galaxy in the Perseus 
cluster JHatch et al J 20061) . This is consistent with an origin for the 
filamentary system being slowly drawn up under the rising radio 
bubbles. 

The spectrum of NGC 4696 exhibits the same strong low ion- 
isation emission lines seen in many other extended emission line 
systems surrounding BCGs. These sy s tems are not all consistent 
with a merger origin. Ijohnstone et al.l j2007h detected strong [Ne 
IIl]A15.56/im emission in their Spitzer IRS SH spectroscopy of 
NGC 4696 and NGC 1275. Opt ical [Ne III] line e mission has also 
been seen ir i both NGC 1275 jHatchet alj|2006l) and weakly in 
NGC 4696 JFarage et allbOlOl) . However, [O III] emission tends 
to be low in the ex t ended emission line sy stems in BCGs (e.g. 
iDonahue et alj200CI) . Ijohnstone et al.l ( 120071) have shown that stel- 
lar photoionisation models struggle to account for the high [Ne 
III] and low [O III] emission unless coupled with high metallic- 
ities. X-ray metallicities of typically solar for oxygen and neon 
dSanders & Fabianll2006h are found in the centre of the Centaurus 
cluster and there is little sign of young stellar populations in the 
filaments. 

The [O IIl]A5007/H/3 emission in NGC 4696 drops dramati- 
cally from the nucleus to ~12 arcsec, indicating the [O III] emis- 
sion is low in the filaments and the majority of emission in the in- 
ner regions is associated with the nucleus. A similar trend is seen in 
many other BCGs (e.g. [jghnstone & Fabian 1988). Our HRB data 
does not extend to beyond the inner filament so we cannot measure 
the [O III] emission in the farthest filaments from the nucleus. The 
ratio of our detected [O III] lines is small compared to shock mod- 
els and the upper limit of the weak [O III]A4364 line indicates a 
[O III] A4364/H/3 ratio of 5 per cent; l ower than the predicted value 
fro m the best velocity shock mo del of lFarage et al.l J201CI) . 

iFerland et alj J2008l l2009h have shown that the characteristic 
LINER spectrum with anomalously low [O III] and high [Ne II] 
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and [Ne III] emission, compared to photoionisation models, can be 
reproduced with their particle heating model. With the exception 
of [N II]A6583, which is under predicted, all our optical emission 
lines detected in the outer filaments in NCG 4696 are within a fac- 
tor of2_or2£ssof their p redicted spectrum. The predicted spectrum 
of iFerlandet alj | |2009|) . uses a gas-phase abundance of nitrogen 
relative to hydrogen of 6.9xl0~^ which correspo n ds to ~0.5Zp) 
using abundances from I Anders & Grevessd l ll989l) . ISanders et al.l 
( 120081) find the nitrogen abundance in Centaurus to be ~4 times 
solar in the central regions, a factor of 8 times that used in the 
predicted particle heating spectrum. The predicted emission line 
ratios do not in general behave linearly with changes in metallic- 
ity of the gas due to the thermostat effect, where increasing the 
abunda nce of a coolant lowers the k inetic temperature to conserve 
energy JOsterbrock & Ferlanal2006i) . However, the emissivities of 
[N II] and Ha are predicted to be relatively constant over a broad 
range of temperatures (10^'^ — 10^ "^ K) within the thermally stable 
regime (see Fig. 17 of iFerland et al.ll200 9) so an increase in the ni- 
trogen abundance by 8 would produce an increase in the [N II]/Hq 
ratio. Note though that this is not necessarily the case for the [N 
I] emission which does not exhibit such a flat emissivity profile. 
Future papers will explore the effect of varying metallicity on the 
particle heating model ratios. 
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APPENDIX A: DATA REDUCTION PLOTS 

Al Sky absorption and emission corrections 

The VIMOS IFU does not have dedicated sky fibres and unlike slit 
spectroscopy an individual spectrum on a fibre does not have both 
regions with pure sky emission and regions with object and sky 
emission. Different fibre spectra are described by different profile 
shape parameters such as the FWHM and skewness. In the wave- 
length direction this results in the the spectral lines having different 
shapes and thus affects the quality of the sky subtraction. An aver- 
age sky spectrum formed by combining spectra with different line 
profiles can result in the presence of s-shape residuals when sub- 
tracted from the data. VIPGI sky subtraction attempts to correct for 
these line profiles by grouping the spectra according to a user de- 
fined sky line. An average sky for each spectrum grouping is then 
calculated and subtracted. The groupings are done on a statistical 
basis thus this sky subtraction method employed by VIPGI is op- 
timised for deep survey observations where the field is devoid of 
extended objects. It is not ideal for observing large galaxies which 
dominate the field of view. For this reason we performed the sky 
subtraction using specific fibres shown to be lacking in emission 
lines associated with NGC 4696. 

Telluric absorption feature corrections for the O2 and H2O 
absorption in the 6000 — 7000 A band was determined from 4 ob- 
servations of standard stars. Standard stars were observed at the 
beginning and the end of the night. The stars chosen are CD-32- 
9927, LTT-7379, Hiltner-600 and LTT-2415. The standard spectra 
G{\) were reduced using VIPGI, as with the science observations. 
These spectra were divided by the ESO spectra Go (A) of the re- 
spective standards, corrected for atmospheric absorption. These ob- 



Figure A2. Top: The contour binning algorithm of JSandersI 120061) ap- 
plied to a 2D brightness map of [N II]A6583 emission in NGC 4696. 
Bottom: The same using the Veronoi tessellation binning algorithm of 
ICappellari & CopinI 120031) . Both maps show the data binned to a signal- 
to-noise of 10. 



servations were then averaged for each VIMOS IFU quadrant and 
the area around the absorption lines were isolated. Fig. lA II shows 
the quadrant one spectra of each of the four standard stars and the 
averaged normalised spectrum. 

Ignoring the curvature of the Earth we recover the object spec- 
trum as seen without the atmospheric absorption (Fo), 



Fo(A) = F(A) 



Go(A) 



G(A) 



(Al) 



for each 15 minute exposure separately so as to correct for airmass 
differences between the standard star (^standard) and the object 

(^object)- 



A2 Propagation of uncertainties 

Our HRO and HRB spectra overlap between 5000 - 6000 A. 
In order to incorporate all our data we split the cubes into three 
wavelength regions. The regions chosen are 4000 — 5100 A, 
5100 - 5800 A and 5800 - 7500 A. The two HR grisms used 
have a slightly different spectral resolution (see Table [TJ with the 
HRO being marginally larger than that of the HRB grism. In or- 
der to amalgamate both data sets it is necessary to re-bin the HRO 
spectra. If we are re-binning original bins i into bins j where some 
fraction Wij of i falls into the jth bin we determine the flux Fj and 
error Sj in the jth bin by, 
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Figure Al. Telluric absoq^tion correction for Ql using four standard star observations. We only correct for atmospheric absorption between 5500 A and 
7000 A. The absorption seen at ~ 5900 A is likely NaD absorption from our Galaxy and as such has not been connected for. The red line shows the median 
observed spectrum of the four standard stars divided by their intrinsic spectrum. The black line is the teluric correction we apply to our data. 



Fj ^Sj ^Yl ^y /' ± 4 / 5Z '^i^^^' 



(A2) 



where fi is the flux in original bin i and ef is the variance in 
bin i. For pixel to pixel fluctuations the error Ej is given by Ej = 

A/X^i '^ii^'i- The HRO and HRB spectra are then summed and the 
total flux Tj in the re-binned grid is given by (Bob Carswell, private 
communicatiorlfl), 



E 



T, 



k{Sj) 



7T2 



E 



(A3) 



kJsW 



where k is the number of spectra summed together. The esti- 
mate of the fluctuations in the spectra are then given by, 



Efe( 



rT2> 



(s^) 



(E 



kjsW) 



(A4) 



this forms the error spectrum used when fitting gaussian mod- 
els to the emission lines and minimising the resulting x^- 



A3 Spatial binning 

To improve the signal to noise in the spectra w e investi g ate th e 
use of two spatial binning algorithms; that of I Sanders! j2006r) . 
which bins the da t a bas ed on surface brightness and that of 
ICappellari & CopinI ( 120031) which uses Voronoi tessellations to pro- 
vide the most 'round' bins. The two techniques applied to our NGC 
4696 data are shown in Fig. lA2l 

Both techniques give similar fits over comparable regions giv- 
ing us confidence that the physical properties of the regions we are 
binning are similar to the limit of our spatial resolution (see Fig. 
|A3l l. Due to the filamentary structure of the optical line nebulos- 
ity surrounding NGC 4696 we choose to bin our spectra using the 
contour binning technique of I Sanders! 1200g). 



ftp://ftp.ast. cam. ac.uk/pub/rfc/vpfit9. 5.pdf 




Figure AS. An example single gaussian fit to the strongest emission lines in 
our spectra. These are the lines of [N II]A6548, Ha A6563, [N II]A6583, [S 
II]A6717 and [S II]A6730. These lines are fit between 6590 and 6820 A. A 
fixed continuum has been subtracted and is estimated as the median value 
between 6700—6750 A. The continuum regions around the Hq+[N II] dou- 
blet and the [S 11] doublet are fit separately. The data and error are shown 
as a thick grey line, the red line indicates the model fit to the data and the fit 
residuals are shown in the bottom panel. 



A4 Velocity component fitting 

The data are fit with multiple velocity components. An example 
single velocity component fit to the Ha, [N II] and [S II] emission 
lines are shown in Fig. |A5l An F-test is then performed to determine 
whether extra velocity components are required. 

The lower panels in Fig. lA4l show the results of the F-test when 
the data is binned by surface brightness contour binning, Voronoi 
tessellations and no-binning. The darker shade implies two veloc- 
ity components are necessary, the light shade indicates the spectra 
were sufficiently well represented with a single gaussian velocity 
component. 



APPENDIX B: EIVIISSION LINE IVIAPS 

The reddening map for a one component velocity fit, on a per pixel 
basis, in the central regions of NGC 4696 is shown in Fig. lBll Con- 
tours of E(B-V) derived from HST B and I band ima ges and Chan- 
dra X-ray Nh column density l lCrawford et al.ll2005n are overlaid. 
IVlaps of our em ission line ratio s, binned using the contour 
binning algorithm of ISandersI yOOg) are shown in Fig. IB2I The 
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Figure A3. FWHM velocity dispersion, in units of kms"^ , for different bi nning teclinique s for a single gaussian fit to the [N 11], Hq and [S II] lines. 
Left: Binning to a signal- to-noise of 10 using the con tour binning technique of lSandersI j2006l) . Middle: Binning to a signal-to-noise of 10 using the Veronoi 
tessellation technique of ICappellaii & CopinI 120031) . Right: The velocity dispersion on a per-pixel basis across the field of view. Here the (0,0) position 
corresponds to the centre of the image. 
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Figure A4. Regions where the fibres are statistically found to be described by one or two velocity components for the specfi'a binned by surface brightness 
(left), Voronoi binning (middle) and no-binning (right). The darker shade indicates two components were necessary to fit the spectra. Regions in a lighter shade 
were found to be modelled sufficiently with a single velocity component. 



top panel showing the single gaussian fits, the middle showing the 
broad central velocity component and the bottom panel in each case 
shows the narrow component. Only pixels where the line emission 
is detected at greater than 4cr are shown. The [O IIl]]/H/3 ratios 
cannot be traced across the whole length of the outer filaments as 
our HRB data only cover the central 27 arcseconds. 



© 0000 RAS, MNRAS 000, 000-000 



20 R.E.A. Canning et al. 



en 

"D 

c 
o 
o 

03 

o 



0.0 



0.4 



0.8 



1.2 



1.7 



2.1 



2.5 




5 

Arc Seconds 



5 0-5 

Arc Seconds 



Figure Bl. Av intrinsic extinction maps derived f rom tlie Ha/H/3 ratio, fit with a single gaussian. Left: Overlaid with contours of E(B-V) Right: Over- 
laid with contours of X-ray Nh ■ Contours are from ICrawford et alj 120051) . The colour bar shows the intrinsic Ay extinction calculated assuming case B 
iOsterbrock & Ferlandl200d) ; only pixels where both Ha and H/3 were detected at 4 sigma or above are shown, the original fits were smoothed by a gaussian 
with FWHM 3 pixels (2"). 
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Figure B2. Clockwise from top left: [N II]A6583/Hq A6563, [S II]A6716/[S II]A6731. [O III]A5007/H/3 A4861, [S II]{A6716 + A6731)/Ha A6563 and [O 
I]A6300/Ha A6563 emission. The top panel in each figure is the ratio in the outer filaments where only one velocity component was detected. The middle 
panel is the narrow component of the two component gaussian fit and the bottom panel is the broad component. 



© 0000 RAS, MNRAS 000, 000-000 



